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Abstract — Network coding and cooperative communication 
have received considerable attention from the research com- 
munity recently in order to mitigate the adverse effects of 
fading in wireless transmissions and at the same time to achieve 
high throughput and better spectral efficiency. In this work, 
we analyze a network coding scheme for a cooperative com- 
munication setup with multiple sources and destinations. The 
proposed protocol achieves the full diversity order at the expense 
of a slightly reduced multiplexing rate compared to existing 
schemes in the literature. We show that our scheme outperforms 
conventional cooperation in terms of the diversity-multiplexing 
tradeoff. 

Index Terms — Cooperative communication, network cod- 
ing, random network coding, outage probability, diversity- 
multiplexing tradeoff. 



I. Introduction 

Channel fading is one significant cause of performance 
degradation in wireless networks. In order to combat fading, 
diversity techniques that operate in time, frequency or space 
are commonly employed. The basic idea is to send the signals 
that carry same information through different paths, allowing 
the receiver to obtain multiple independently faded replicas of 
the data symbols. Cooperative diversity tries to exploit spatial 
diversity using a collection of distributed antennas belonging 
to different terminals, hence creating a virtual array rather than 
using physical arrays. 

In iXi Ahlswede et al. introduced network coding to achieve 
the max-flow value for single-source multicast which could 
be impossible by simply routing the data. Since then, network 
coding has been recognized as a useful technique in increasing 
the throughput of a wired/wireless network. The basic idea of 
network coding is that an intermediate node does not simply 
route the information but instead combines several input pack- 
ets from its neighbors with its own packets and then forwards 
it to the next hop. However, since network coding is devised at 
the network layer, error-free communication from the physical 
and medium-access layer is usually assumed, which is a 
simplifying assumption for wireless communications. 

Efforts have also been made to apply network coding to 
the physical layer, e.g. in |9|. Towards that goal, cooperative 
schemes have been proposed that make use of network coding 
in a cooperative communication setup, and studies have been 
conducted to determine whether network coding provides any 
advantages over existing cooperative communication tech- 
niques ii, im. 




Fig. 1. System model: A' source-destination pairs and M relays 



In J&l, a network-coded cooperation (NCC) was pro- 
posed and its performance was quantified using the diversity- 
multiplexing tradeoff analysis which was originally proposed 
for multiple antenna systems in ifTol . NCC was shown to 
outperform conventional cooperation (CC) schemes which 
includes space-time coded protocols |4| and selection relay- 
ing 10 • NCC requires less bandwidth, and yields similar 
or reduced system outage probability while achieving the 
same diversity order However, these results are based on an 
optimistic assumption that any destination node should receive 
the packets that are not intended for it without any error so 
that the intended packet can be recovered from the xor'ed 
packet sent by the relay. When this assumption is removed 
the scheme can no longer achieve the full diversity order of 
M + 1, where M is the number of cooperating relays, but only 
a reduced diversity order of 2. 

In this paper, we propose a network coded cooperation 
schemes for N source-destination pairs assisted with M relays. 
The proposed scheme allows the relays to apply network 
coding on the data it has received from its neighbors using 
the coefficients from the parity-check matrix of a {N + 
M,M,N + 1) MDS code. A closed form expression for the 
outage probability is derived. We also obtain the diversity- 
multiplexing tradeoff performance of the proposed scheme. 
Specifically, it achieves the full-diversity order M + 1 at the 
expense of a slightly reduced multiplexing rate. We show 
that our scheme outperforms NCC and CC in terms of the 
probability of outage. 

The rest of the paper is organized as follows. Section 
Hn discusses the system model, description of the proposed 
scheme and the design of the network code. In Section |llll 
we present the main result and the proof. In Section |IV] the 
performance of the proposed scheme is compared in terms of 
DMT and average outage probability. Section [V] contains the 
conclusions. 
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Fig. 2. time-division allocation for the different schemes compared: (a) CC 
(b) NCC (c) DNCC 



II. System Model 
A. General System Description 

The network studied in the paper is composed of N source- 
destination pairs denoted as (si, di), . . . , (sat, dAr), and M 
relays denoted as ri , . . . , tm in a single-cell where all the 
nodes can hear the transmissions of each other as shown in 
Fig. [T] We consider two different transmission scenarios. In 
the first scenario, each source node s, is trying to transmit the 
data packet 6i to all the destinations di ,i ^ 1, . . . , N which is 
known as the multicast scenario. In the second scenario, each 
source node Si is trying to transmit the data packet di to only 
destination di and we will refer to this scenario as the unicast 
scenario. We assume that each packet is composed of L bits: 

= 9i^2, ■ ■ ■ , ^i.i]- As in |3|, these length-i blocks of 
bits that are transmitted on each link will be treated as elements 
of a finite field F^, g = 2^. All the nodes are assumed to be 
equipped with half-duplex (i.e. cannot transmit and receive at 
the same frequency) single-antennas. Each data packet Qi is 
error control coded and modulated, and transmitted in T time 
slots. 

We assume the users transmit in equal time slots with 
transmit power P. The channel between any pair of nodes is 
assumed to be frequency flat fading with additive white Gaus- 
sian noise (AWGN). Let Xj e C^^"^ denote the transmitted 
symbols from node j and S C^^"^ the received symbols at 
node i, and the additive noise has independent and identically 
distributed (i.i.d.) entries z^.f ^ CN(0, 1), and hij E C the 
instantaneous channel realization. Then, the channel within 
one block can be written as 



(1) 



where p is the average received SNR at the destination. In the 
above equation, the transmitter could be any of the sources or 
relays, the receiver could be any of the relays or destinations, 
as long as the transmitter and receiver are different (i.e., not 
the same relay). We assume that the channel coefficient hij 
remains constant during the transmission time of a packet. 

The channel coefficient hij between any two nodes is 
modeled as i.i.d. with zero-mean, circularly symmetric com- 
plex Gaussian random variables with common variance 1/(3. 



Therefore, |hi are exponentially distributed with parameter 
13. A total of NL bits are transmitted in {N + M)T channel 
uses, therefore the system rate is i? = (n+m)t ^^'■^ P^'" 
channel use (BPCU). The transmission rate Ri for one source 
or one relay per one packet is fixed, identical, and equal to 
Ri^k^ E±^R BPCU. 



B. Deterministic Network Coded Cooperation (DNCC) 

Our transmission scheme consists of two stages; see 
Fig. |2|c). In the first stage, direct transmissions from the 
sources to the destinations take place in N orthogonal time 
slots. Thanks to the broadcast nature of the wireless medium, 
all the destinations and the relays overhear the transmissions. 
At the end of the first stage, each relay tries to decode all 
N packets. If a relay can successfully decode all the packets 
(we denote this assumption by si/), then it participates in the 
second stage. Otherwise, it remains silent. 

In the second stage, the participating relays perform network 
coding. Specifically, relay i will transmit the linear combina- 
tion 
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(2) 



where Qk is the corresponding finite field element in Fg for 

Ok- 

The coefficients a^'s are predetermined and they are de- 
signed in a way to maximize the chance that the received 
linear combinations are actually decodable at the destination. 
We discuss the problem of how to choose these predetermined 
coefficients in detail in Sec. III-CI 

In order to express the overall transmitted data packets, we 
define the following matrix: 
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(3) 



where the top part is an iV x identity matrix. We also define 
the iV X 1 finite field vector corresponding to the original 
source packets as 



e = [91,62,. 



(4) 



where {■)'^ denotes transpose. Using matrices A and Q, the 
potential transmitted data packets can be expressed as X — 
A&. 

C. Design of the Linear Network Coding Matrix 

Each destination di will attempt to decode all the packets 
that have been transmitted at the two stages. It is possible 
that some of the packets cannot be decoded due to severe 
channel fading. However, if enough number of packets can be 
decoded, and the linear network coding coefficients a^.n are 
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Fig. 3. DMT comparison of various schemes 

judiciously designed, then the destination will be able to obtain 
the transmitted packets by solving a linear system of equations. 
In this subsection, we discuss the issue of the design of the 
coding matrix A (or the bottom M rows of A to be exact). 

The row Kruskal-rank Q of A, denoted as k{A), is the 
number r such that every set of r rows of A is linearly 
independent, but there exist one set of r + 1 rows that are 
linearly dependent. The column Kruskal-rank can be defined 
similarly. 

The following result relates the column Kruskal-rank of the 
parity-check matrix of a linear block code to its minimum 
distance dmin'- The column Kruskal-rank k{H) of the parity 
check matrix H of a linear block error-control code and the 
minimum Hamming distance dmin of the code are related by 
dmin — i^{H) + 1. Please refer to fSf for a proof. 

We would like to have the matrix A to have a large row 
Kruskal rank so that a minimum number of equations is 
needed by any destination to solve for the source packets. 
The row Kruskal rank is of course less than or equal to the 
rank, which is N. To maximize the Kruskal rank, we should 
design A such that it has Kruskal rank equal to N. Or stated 
differently, any N rows of the matrix A should be linearly 
independent. Depending on the sizes N and M, this may or 
may not be possible in a given finite field F^. Note that whether 
this is possible is closely related whether maximum distance 
separable (MDS) codes exists in a certain field for certain code 
dimensions [5]. 

The transpose of the parity check matrix of a systematic 
{N + M, M, + 1) MDS code can be used as an encoding 
matrix A for our DNCC scheme to minimize the total number 
of packets necessary at the destinations for decoding the source 
packets. If such an A is used, then each destination needs and 
only needs packets (from the sources and relays) for correct 
decoding. 

III. Diversity-Multiplexing Tradeoff 

As mentioned in the introduction, we will investigate 
the performance of the proposed scheme via diversity- 
multiplexing tradeoff (DMT). DMT is widely accepted as a 
useful performance analysis tool in cooperative systems [2j, 
||6l. For completeness, we give the formal definitions as in 
fW\: A scheme C{p) is said to achieve spatial multiplexing 
gain r and diversity gain d if the data rate is 

Rip) 



p->oo log(p) 



(5) 



and the average error probability is 

l og(^e(p)) 

log(p) 



lim 



= -d. 



(6) 



A. Outage Probability 

In ifTO l, for long enough block length T it was shown that 
the probability of error is dominated by the outage probability. 
Therefore, we will only consider outage probabilities in our 
analysis. The instantaneous mutual information of the channel 
model in ([T]i is given by: 



/(X,;Yj) -log(l + |h,,,fp). 



(7) 



where and Yj denote the transmitted symbol by node i 
and received symbol by node j. We write a(r) = 6(r) if 
linir^oiO'i'T) /b{T)] = 1. An outage event between a trans- 
mitter i and receiver j occurs when the instantaneous mutual 
information for the channel hij is less than the transmission 



rate Ri BPCU. Since |hi ., | is exponentially distributed, the 
outage probability for the channel in ([T]) is given by: 



Po=Pr{I{X^:Y/}<R,) 



Pr(|h,,,f < r) = 1 - exp(-/3r) ^ /3r, 



(8) 



where r = 



B. Further Improvements 

1) Decoding at the relays: Decode-and-forward schemes 
suffer from performance loss when the source-relay channel 
is in outage. And if a multi-source scenario is considered the 
performance loss becomes even more severe. Therefore, the 
assumption that the relay has to decode all the packets in 
order to be able to cooperate becomes a bottleneck for such 
schemes. We could relax the assumption jz/ and assume that 
the relays will participate cooperation even though they have 
not been able to decode all the packets. Specifically, in the 
second stage, the participating relays perform network coding 
on the packets that they have received correctly. If relay i was 
able to decode the packets correctly from the sources in the 
set Si where Si C {!,..., A^}, then it will transmit the linear 
combination 

a^^kOk (9) 

keSi 

where 9a; is the corresponding finite field element in for 

9k- 

Theorem 1. The diversity-multiplexing tradeoff of the linear 
network coded cooperation with M intermediate relay nodes 
for both multicast and unicast and both either with or without 
the assumption is: 



d{r) = (M + l) 



1 



(N + M) 

1 

N 



,re 0, 



N 



N + M 



(10) 



Proof: Due to severe fading some of the received packets 
may not be successfully decoded by the destination. This 
can be viewed as some of the rows of A being erased. We 
denote the resulting submatrix by Ai for destination di. The 
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successfully received packets at destination di after decoding 
can be expressed as Yi = AiQ. In the multicast problem, 
destination di cannot recover 8i when the submatrix Ai is 
rank deficient, i.e. when rank{Ai) < N. This happens when 
less than channels are not in outage which results in an Ai 
matrix that has at most iV — 1 rows. We define this event as 
Eli — {Ai has at most — 1 rows}. 

Next, we define several other events which will simplify the 
analysis. Let Em denote the event that m relays fail to receive 
all the 8i's correctly. Then we have: 



we have: 



N 



P{£) = Y[Pr{IsAX;Y)> R,) (11) 

i=l 
N 

= l[Pr{\h,A^ >t)=cxp{-N(3t) (12) 



Similarly, define E{m,n) to be the event that n channels 
out of ni were in outage: 



P{E{m,n)) 



Poil-PoY 



(13) 



where Pq is given by (|8]l. Now, we can express P^Eu) as: 

M N+M-m 

P{Eu)=Y.P(^ra)- P{E{N + M-m,n)) 

m— n— M— m+1 

(14) 

The first summation stands for the probabiUty of the event 
that m of the relays were in outage, leaving us with only 
N + M — m nodes which were not in outage for that particular 
destination. Notice that, as p — > oo, r — > 0. We would like to 
find the following limit: 



lim 



P{Eu) 



(15) 



We consider the individual terms in the summations one-by- 
one and find the term with the smallest order of r. Observe 
that P{Em) = KmT"^ where Km is a constant: 



• lim 



(1 - exp{-Nl3T)y 



M 



Next, consider the term P{E{m,n)), 



r-i-O t" 



m 
n 



Po"(l - Po)'"~" 



/3" 



(17) 



When n is equal to M — m + 1, we have the smallest order 
r term . Substituting n = M-m+1 in P{E{N + M -m,n)). 



lim{ ( ^ + ~ ? ) (1 - p,)^+M-m-iM-m+l) 

T^o \ M—m + 1 ' 

P^'-"'+\ N + M^m\ M-m+i 

_A/-m+i / [ M -m+1 r' 



Therefore, we have: 

^. PiEm) P{E{N + M-m,n)) 

'"Vivffl-r";"-")^"— (18) 

\ Ai — TO + 1 / 

Finally, using ( fTsT i, we have: 



M -rn+1 



'N + M_„_, \ M+1 J(J 
Z N \ 



P 



(19) 



Now, suppose the fixed rate has chosen to be: 

i? = r log p (20) 
Substituting ( l20l i into ( fT9] l, we arrive at: 

-/MWiV + Af-TOX 

Since the summation term is just a constant in terms of p, we 
have the desired result. We skip the proof for the unicast and 
without assumption cases due to space limitations. ■ 

IV. Comparison WITH Other Schemes 
A. DMT Comparison 

In this section, we would like to compare the proposed 
scheme with the existing schemes in the literature in terms 
of diversity-multiplexing tradeoff. The closest scheme in the 
literature is the Network-coded cooperation (NCC) considered 
in |6|. There are N s-d pairs and M relays in a single cell. 
However, instead of all the relays, only one relay transmits 
following the direct transmissions from the source nodes, 
which results in total of + 1 time slots (Fig. |2] (b)). Using 
fewer time-slots NCC achieves a better spectral efficiency than 
DNCC. However, NCC can only provide a fixed diversity order 
of 2, while the proposed scheme achieves the full-diversity 
order of Af + 1. 

In the following, for comparison we include the DMT 
performance of NCC and that of conventional cooperation 
(CC). 

The diversity-multiplexing tradeoff of NCC is given by f6l: 



d(r) = 2(l-^r 



e 0, 



N 



N + 1 



(22) 



The DMT of the CC schemes with M intermediate relay 
nodes is given by E): 



d(r) = (Af + 1) (1 - 2r) , rS (0,0.5). 



(23) 
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We present diversity-multiplexing tradeoff of the existing 
schemes and DNCC in Fig. [3] As can be seen from the figure, 
although DNCC can achieve the full diversity order, NCC fails 
to do so. We also see that DNCC achieves a better DMT than 
CC when N > M. 

B. System Outage Probability 

Here, we compare the system outage probability of DNCC 
with the other schemes. The system outage occurs when any 
di is unable to decode 8^ reliably: 

N 

Ps = l P{Eu)) (24) 

1=1 

We compare (l24l i with the system outage probabilities (30), 
(43) derived in |6|. In Fig. |4(a)| the proposed method clearly 
outperforms NCC by achieving the full diversity of M + 1 as 
compared to NCC's fixed diversity order of 2. 

C. Monte-Carlo Simulation 

We also compare the proposed scheme with the existing 
schemes via Monte-Carlo simulations. In the simulations, 
only channel conditions are considered to isolate the diversity 
benefits of the scheme. We generate an {N + M) x N and an 
N X N matrix that contains the channel coefficients for each 
destination and each relay, respectively. Then, we decide that 
the transmission is successful for any link if the instantaneous 
channel condition is large enough to be able to support the 
given data rate and we update the same size linear coefficient 
matrices accordingly. After all the transmissions take place, we 
perform Gaussian elimination on the updated linear coefficient 
matrices to conclude whether each destination di was able to 
recover the source packet Qi or not. The channel coefficient 
variances are chosen to be equal to one and Ri — 1 BPCU. 
Please note that we considered the average outage error 
probability which is found by dividing the total number of 
errors occurred by the number of source nodes instead of the 
system error probability. In all the figures only the unicast 
scenario is adapted since CC cannot be implemented in a 
multicast scenario. As it is proven by the Theorem [T] the 
performance loss incurred due to the assumption jz/ is not in 
terms of diversity gain but it is in terms of coding gain. This 



is validated through simulations as shown in the Fig. |4(b)| and 

Fig. EMI 

V. Conclusions 

In this paper, we have proposed a network coded coop- 
eration scheme for N source-destination pairs assisted with 
M relays. We studied two different transmission scenarios, 
unicast and multicast. The proposed scheme allows the relays 
to apply network coding on the data it has received from 
its neighbors. We establish the link between the parity-check 
matrix for a {N+M, M, N+1) MDS code and the coefficients 
to perform network coding in a cooperative communication 
scenario consisting of N source-destination pairs and M 
relays. We also obtained the diversity-multiplexing tradeoff 
performance of the proposed scheme, and showed its advan- 
tage over the existing schemes. Specifically, it achieves the 
full-diversity order M + 1 at the expense of a slightly reduced 
multiplexing rate. 
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